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I .  INTRODUCTION 


Monte  Carlo  dose  attenuation  kernels  for  the  transport  of 
fission  electrons  through  aluminum  have  been  calculated  for 
simple  slab  and  spherical  geometries.  These  kernels  are  pre¬ 
sented  in  Section  II. 

Solid  angle  sectoring  methods  are  commonly  used  for  complex 
geometry  calculations.  The  selection  of  kernels  appropriate  to 
sectoring  calculations  is  discussed  in  Section  III. 

Several  factors  of  2  are  sometimes  handled  incorrectly 
in  space  radiation  dose  calculations.  These  factors  and  some 
mathematical  relationships  are  presented  in  Section  IV. 

Section  V  derives  several  relationships  for  spherical  shields 
and  describes  sectoring  kernels  that  predict  correct  results  in 
two  simple  shield  geometies. 

Section  VI  reiterates  prior  information  on  sectoring  kernels 
and  includes  a  scheme  for  predicting  shield  curvature  effects. 

Dose  attenuation  data  for  slab,  sphe?  •  '  and  solid 

sphere  shields  are  presented  in  the  appendix.  These  data  span 
the  energy  range  of  0.1  to  8  ileV  and  materials  from  hydrogen 
through  uranium. 
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II.  KERNELS 


Attenuation  kernels  are  denoted  as  functions  of  shield 

thickness  t;  e.g.,  K(t)  is  the  response  behind  a  shield  of 

thickness  t.  The  units  of  t  are  given  variously  as  cm,  mils, 
2 

and  g/cm  , 

All  kernels  are  presented  in  the  form  used  for  sectoring 
analyses.  If  R  denotes  the  total  response  e.g.,  dose 


R 

complex 

geometry 


solid  4* 
angle 
sectors 


K(t  of  sector) 


(1) 


This  implies  that  slab  shield  kernels  which  have  a  source  inci¬ 
dent  on  only  one  side  are  multiplied  by  2  to  form  K(t).  Thus, 
source-shield  configurations  have  mirror  symmetry  in  the  dose 
plane. 

For  example,  a  dose  point  between  two  slab  shields  of 
thicknessess  t^  and  tg  would  have  the  response  computed  as 

Rtotal  “  *  +  <2> 

and  if  t,  becomes  infinite 

“total  *  »  «‘l>  <3> 


The  geometry  of  the  shield  is  denoted  by  a  subscript  on 
the  kernel  X: 

K_(t)  is  two  times  the  dose  behind  a  slab  shield  of  thick¬ 
ness,  t.  No  backscatter  material  is  present.  Radiation  is 
incident  from  one  side  only. 

X_(t )  is  two  times  the  dose  behind  a  slab  shield  of  thick- 
00 

ness,  t,  with  an  infinite  backscatter  material.  Radiation 
is  incident  from  one  side  only. 

K_(t)  is  the  dose  between  two  slab  shields  each  of  thick¬ 
ness,  t.  Radiation  is  incident  from  both  sides.  The 
symmetry  makes  Kja(t)  equivalent  to  two  times  the  dose  behind 
a  shield  of  thickness,  t,  with  a  backscatterer  of  thickness 
t. 

XQ(t)  is  the  dose  at  the  center  of  a  spherical  shell  shield 
of  thickness,  t.  The  inner  radius  of  the  spherical  shell  is 
essentially  infinite.  Radiation  is  incident  on  the  outside 
of  the  shell.  Radiation  transmitted  through  the  shell  can 
baoiscatter  from  the  shell. 

K.(t)  is  the  dose  at  the  center  of  a  solid  spherical  shield 
of  thickness,  t,  with  radiation  incident  on  the  outer  surface 
of  the  sphere. 

The  dose  for  zero  shield  thickness  is  the  same  for  geometries 
without  a  backscatterer 

KJO)  -  K.<0)  -  ko<0)  -  K. (0)  -  Ktrae  splce  (4) 
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The  double  slab  shield  configuration  is  equivalent  to  an  infinite 
backscatterer  for  large  thickness  so  that 

K_(t)  »  K,(t)  for  large  t  (5) 

00 

CALCULATION  TECHNIQUES 

An  isotropic  radiation  field  0°  in  particles/cm2  is  described 
as  a  transport  problem  by  a  surface  source  or  current  of  strength 

.  n 

jr/4  particles/cm.  The  surface  on  which  the  source  is  defined 
must  be  exterior  to  the  shield  geometry.  A  cosine  angular  distri¬ 
bution  relative  to  the  inward  surface  normal  is  required  to 
reproduce  an  isotropic  flux  inside  the  source  surface. 

A  fission  electron  energy  distribution  was  used  in  the  dose 
calculations.  The  energy  variation  of  this  spectrum  was  repre¬ 
sented  by  (normalizing  constants  suppressed) 

f (E)  «  exp(-  0.575  E  -  0.055  E2)  particles/MeV 


where  E  is  the  electron  energy  in  MeV,  restricted  to  the  range 

0  <  E  <  7  lieV,  This  energy  dependence  was  integrated  numerically 

2 

for  normalization  to  a  fluence  of  one  electron/cm  . 

The  Monte  Carlo  calculations  were  preformed  using  adjoint 
methods  described  in  Reference  1.  Histories  were  generated  for 

T.  Jordan,  TM. ,  An  Adjoint  Charged  Particle  Transport  Method. 
EMP.L7n.072,  July  l9t6. 


an  infinite  medium  of  aluminum.  Shield  geometries  of  various 
thicknesses  were  o/erlaid  onto  the  infinite  medium  electron 
tracks  to  obtain  a  dose  contribution  for  each  geometry  (Ref.  2). 

The  production  calculation  required  5  minutes  (CDC  7600)  to  obtain 
simultaneous  results  for  the  five  geometries  and  40  shield  thick- 
nessess.  Calculated  standard  deviations  less  than  2  percent 
were  obtained  for  most  data  points. 

CALCULATION  RESULTS 

Dose  rates  are  expressed  as  point  results  in  rad  (Si)  units. 

No  silicon  was  actually  included  in  the  ilations.  The  flux 
spectrum  obtained  in  aluminum  was  LET  weighted,  i.e..  folded 
with  the  restricted  linear  energy  transfer  function  for  silicon. 

The  restricted  energy  transfer  excluded  energy  loss  mechanisms 
included  in  the  transport  simulation,  i.e.,  radiative  energy  loss 
and  energy  losses  corresponding  to  high  energy  secondary  electrons. 

The  effect  of  including  a  backscatter  material  in  dose 
calculations  is  shown  in  Figure  3 .  The  lower  curve  is  for  a 
one  side  slab  without  a  back^catterer .  K_(t),  and  the  upper 
curve  is  for  the  two  slab  shield  K^t),  i.e.,  backscatter 
thickness  equal  to  shield  thickness.  For  shield  thicknesses 
greater  than  10  mils  aluminum,  the  b&ckscattered  electrons  add 
about  50  percent  to  the  dose.  Restated,  the  absence  of  a 
backscatte.’  material  reduces  the*  dose  by  approximately  35  per¬ 
cent  . 

2.  Jordan,  T.SJ. ,  Parametric  Adjoint  Electron  Transport  in  1-P 
Geometries,  KIP. L76.G97 .  September  l 97*!l 


Tbe  effect  of  infinte  backing  is  mnst  notable  at  zero 
thickness  where  the  slab  dose  is  increased  by  approximately 
30  percent.  This  is  indicated  by  tbe  highest  curve  in  the 
vicinity  of  zero  shield  thickness  The  infinite  backing  curve 
cannot  be  distinguished  from  the  double  shield  curve  for  thick-* 
nesses  greater  than  25  oils. 

The  effect  of  shield  geometry  on  dose  is  shown  in  Figure  2- 
The  lowest  curves  are  for  the  slab  geometry  shields  which  have 
backing. 

Tbe  intermediate  curve  is  for  spherical  shell  shields. 

The  shell  shield  curve  roughly  parallels  the  slab  shield  curve 
for  thicknesses  of  30  mils  or  more.  The  shell  shield  dose  is 
higher  by  30  to  50  percent.  The  parallel  nature  of  the  slab 
and  s  ell  curves  implies  that  the  angular  distribution  of  the 
electrons  has  an  approximate  equilibrium  since,  for  an  infinite 
spherical  shell  radius 

K0(t)  «*  4*Ku(t,l)  (7) 

where  K„(t,u)  is  the  dose  per  steradiao  for  slab  geometry  as  a 
function  of  u,  tbe  cosine  of  the  angle  between  the  slab  normal 
and  the  particle  direction. 

The  highest  dose  curve  is  for  a  solid  sphere  shield  geometry. 
For  this  geometry,  the  dose  curve  diverges  from  the  other  curves 
as  thr.  shield  thickness  increases. 
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The  enhanced  dose,  or  focusing,  results  from  the  shorter 
average  particle  paths  traversed  in  the  spherical  geometry. 

More  of  the  external  source  electrons  are  able  to  penetrate  the 
shield  in  the  solid  sphere  geometry.  As  the  shield  thickness 
approaches  the  end  of  the  effective  particle  range,  only  a  small 
•spot  -of  aran  or.  Jthe  slab  shield  surface  can  contribute  to  the 
dose.  For  the  same  shield  thickness,  the  entire  surface  of  the 
solid  sphere  will  contribute  to  the  dose. 

ACCURACY 

The  accuracv  of  these  curves  is  determined  partially  by 
comparisons  witn  other  calculations.  Figure  3  is  a  comparison 
of  the  infinite  backed  dose  (slab  shield)  with  Morel’s  Monte 
Carlo  calculations  reported  in  Reference  3. 

The  calculated  results  in  Reference  3  have  been  multiplied 
by  2  in  correspondence  with  the  terminology  used  here  (all  one 
side  slab  results  multiplied  by  2  to  correpond  with  4ir  source 
incidence).  These  data  are  for  the  average  dose  in  a  silicon 
device  1  mm  thick  with  infinite  backing  in  contrast  to  the  point 
dose  results  reported  earlier.  Silicon  and  aluminum  have  very 
similar  cross  sections.  Therefore,  the  1  mm  Si  is  comparable  to 
23  mils  of  aluminum.  This  average  dose  data  has  been  plotted 
in  histcrgram  form. 

3.  Morel,  J. E. ,  Doses  to  a  Thin  Silicon  Slab  Behind  Aluminum, 
Alumlnum-Tantalum,  and  Aluminum-Lead  Shields  for  Isotropic 
Fission  Electrons',  AFWL-bYT-TN-^S-l ,  Air  Force  Weapons  “““ 
Laboratory i  Kirtlani  AFB,  NMex,  July  1975. 
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Figure  3. 
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Fission  Spectrum  Electrons,  Unit  Free  Space  Fluence 


There  is  a  definite  divergence  of  the  calculations  which 
approaches  30  percent  at  175  to  200  mils.  The  reason  for  the 
divergence  cannot  be  identified  from  material  at  hand. 

The  experimental  results  shown  in  Reference  4*  have  been 
used  in  previous  validations  of  the  adjoint  method.  These 
results  were  obtained  from  hollow  boxes  of  various  aluminum 
wall  thicknesses.  The  data  for  30,  60,  and  124  mil  walls  are 
used  for  comparison  by  the  following  steps:  near  and  far  wall 
doses  are  added,  then  'he  sum  is  doubled.  The  experimental 
data  are  for  a  source  incident  on  one  wall.  Adding  near  and 
far  wall  dose  accounts  approximately  for  backscatter.  Doubling 
accounts  lor  two  side  incidence. 

BREUSSTRAHLLNG 

The  infinite  medium  tracking  models  used  for  these  calcula¬ 
tions  did  not  produce  sufficient  photons  to  yield  accurate  photon 
dose  results  in  the  small  shield  g^om^tnes.  Peripheral  numer¬ 
ical  calculation  indicates  a  peak  bremsstrahlung  dose  of  npprox- 
-12 

imately  f  x  1C  rad  (Si)  por  unit  free  space  electron  fluence 
in  the  slab  geometries. 


4.  Van  Gunten,  0.0. ,  Three  Dimensional  Fiectron  Dose  Distribu¬ 
tion  Measurements,  University  of  Maryland  Thesis.  197 i. 

*  Van  Gunten,  0.0.,  Private  Communication  on  reduction  of 
experimental  data  by  two,  February  1976. 
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III.  SECTORING 


i 

i 


1 


Complex  geometry  dose  calculations  can  be  performed  effi¬ 
ciently  by  adjoiut  Monte  Carlo.  The  most  usual  calculation  method, 
however,  is  solid  angle  sectoring. 

Sectoring  calculations  must  neglect  details  of  the  particle 
transport.  However,  these  approximate  methods  should  yield  the 
correct  answer  when  applied  to  simple  geometries. 

The  usual  sectoring  method  is  to  integrate  numerically 

R  =  4i  J  dH  (8) 

4  v 

where  t(u)  is  the  mass  thickness  encountered  between  the  detector 
point  and  the  exterior  of  the  vehicle  along  the  direction  u. 

The  K(t)  is  a  response  attenuation  kernel  either  tabulated  or  curve 
fit  as  a  function  of  mass  thickness  t. 

By  requiring  that  sectoring  methods  yield  the  correct  dose 
at  the  center  of  a  solid  sphere  of  radius  t,  it  follows  that 
K(t )  must  be  the  solid  sphere  kernel  K.(t). 

Requiring  that  sectoring  also  produce  the  correct  dose  at 
the  center  of  a  semi-infinite  slab  of  thickness  2t,  it  follows 
that 

1 

R  *  Ka(t)  «  i  K.<t/y)  dy  (9) 
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Calculating  partial  derivatives  with  respect  to  t  (easier  if  the 
integration  variable  is  first  changed  to  r  =  t/y)  yields 

dK_(t) 

K.(t)  =  K,(t)  -  t  — gt  -  (10) 

This  relationship  must  be  satisfied  by  the  kernels  if  sectoring 
analysis  gives  the  correct  answer  for  the  limiting  cases  of 
simple  slab  or  solid  sphere  shields.  The  relationship  has  been 
reported  in  a  private  communication.* 

The  relationship  is  exact  for  cosine  sources  and  the 
straight  ahead  approximation,  i.e.,  particle  transport  with 
energy  loss  but  no  angular  deflections.  The  straight  ahead 
approximation  is  applicable  to  and  used  for  heavy  charged  parti¬ 
cles.  Application  to  electron  problems  is  suspect  for  energies 
less  than  tens  of  MeV. 

Equation  10  is  also  exact  for  the  following  hypothetical 
transport  problem:  a  uniform  Infinite  medium  and  a  cosine  angular 
distribution  surface  source  which  emits  positive  (real)  particles 
into  the  hemisphere  towards  the  detector  and  negative  (not  real) 
particles  into  the  hemisphere  away  from  the  detector. 

The  hypothetical  problem  has  two  aspects  not  encountered 
in  real  problems.  Source  particles  emitted  into  the  forward 
(detector)  hemisphere  may  make  multiple  passes  across  the  source 
plane,  enhancing  response  kernels.  Negative  particles  emitted 


*  Smith,  E.C.,  Private  Communication,  Hughes,  1976. 
Radke,  G.,  Private  Communication,  AFWL,  1977. 


into  the  backward  hemisphere  may  also  cross  the  source  plane 
one  or  more  times,  contributing  and,  thereby,  decreasing  the 
kernel.  The  net  effect  is  a  kernel  decrease. 

Arguing  that  particles  reflected  across  the  source  plane 
are  a  second  order  effect  gives  credence  to  using  the  relation¬ 
ship  generally.  More  support  is  obtained  by  checking  the 
equation  for  problems  where  both  slab  and  sphere  kernels  are 
available. 

Slab  and  sphere  kernels  have  been  calculated  for  Jovian 
electron  spectra.  Testing  of  Equation  10  for  Jovian  kernels 
is  shown  in  Figure  4.  The  agreement  is  good  for  both  aluminum 
ana  tungsten.  However,  the  Jovian  spectrum  has  appreciable 
content  above  10  MeV  where  the  straight  ahead  approximation 
also  starts,  to  apply. 

Equation  10  can  be  tested  on  the  fission  electron  kernels 
with  the  results  shown  in  Figure  5,  Again,  the  agreement  is 
good.  The  fission  electron  spectrum  is  harder  than  radiation 
belt  spectra.  However,  it  has  sufficient  low  energy  content, 
less  than  2  MeV,  to  suggest  that  Equation  10  holds  under  fairly 
general  conditions. 

The  following  steps  are  suggested  for  sectoring  analysis: 

(1)  evaluate  K_(t);  i.e.,  a  Monte  Carlo  calculation  for 

a 

a  source  emitting  J  particle  with  a  cosine  distribution 
in  the  forward  hemisphere  incident  on  a  slab  shield  of 
infinite  thickness. 
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(2)  calculate  K.(t)  using  Equation  10; 

(3)  use  K.  (t)  in  the  sectoring  analysis. 

The  infinite  slab  thickness  ensures  that  backseattcred 

particles  are  accounted  for  in  the  analysis.  The  kernel  K_(t) 

00 

is  calculated  as  basic  data  since  almost  all  Mcnte  Carlo  codes 
can  compute  it.  Nonadjoint  methods  have  difficulties  when 
computing  K.(t)  directly. 

The  relationship  between  slab  and  sphere  kernels  applies 
to  solid  spheres  and  not  spherical  shells.  From  the  adjoint 
calculations  for  spherical  geometries,  it  is  seen  that  the 
solid  sphere  response  is  greater  than  the  spherical  shell 
response.  However,  the  solid  sphere  kernel  will  yield  solid 
sphere  response  at  the  center  of  a  spherical  shell.  This  over 
estimate  for  one  simple  geometry  is  deemed  better  than  using 
kernels  which  under  estimate  the  response  in  other  simple 
geometries. 

A  more  complicated  sectoring  kernel  is  given  in  Section  V. 
This  kernel  requires  modifications  to  the  sectoring  program. 

The  kernel  eliminates  the  conservatism  of  the  solid  sphere  kernel 
when  applied  to  shell  geometries  and  produces  the  correct  response 
when  applied  to  simple  geometry  solid  shields,  slabs,  and  the 
interior  of  spherical  shells. 


IV.  MANIPULATIONS 


A  space  radiation  environment  is  usually  expressed  as  a 
scalar  fluence  (or  flux)  differential  in  energy 

$(E)  (particles/cm2  MeV)  (11) 

or,  equivalently,  by  the  integral  scalar  fluence 


where 


$(E’)  dE' 


(12) 


♦(E)  *  -  {jg  ( >E) 


(13) 


This  fluence  environment  has  inverse  area  units  which  is 
often  misleading.  The  density  per  unit  area  is  measured  on  a 
surface  perpendicular  to  the  direction  of  particle  motion.  A 
scalar  fluence  implies  integration  over  all  particle  directions. 
Therefore,  the  density  per  unit  area  of  a  surface  with  fixed 
orientation  cannot  always  be  determined. 

With  isotropy,  the  differential  angular  flux  can  be  written 
as 

$(x.u,E)  «  £i|2  particles/cm2  MeV  steradian  second  (14) 

where  x  denotes  position  and  u  is  a  unit  vector  in  the  direction 
of  particle  motion. 
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A  more  basic  definition  of  the  differential  angular  flux 
is  the  particle  density 


0(x,u,E) 


(particles/cm  MeV  steradian)  (15) 


multiplied  by  the  particle  speed  v(E)  in  cm/s 


<J>(*,u,E)  •  o(x,u,E)  v<e> 


(16) 


The  vector  current  is  defined  as  the  flux  times  the  particle 
direction  vector  u 


i(x.u.E)  *  u  $(x.u,E) 


(17) 


the  differential  vector  current  can  also  be  written  as  the 
particle  density  times  the  vector  velocity  v(E)  »  u  v(E) 


,J(x.u,E)  «  v(E)  p(x.u.E) 


(16) 


Defining  a  fixed  unit  area  in  space  and  letting  n  denote 
the  unit  normal  to  this  area  (on  the  side  containing  the  pro* 
jection  of  u),  then  the  number  of  particles  crossing  this  unit 


area  is  the  component  of  the  current  ^  along  n 

*  H-,l(s.u.£) 

-  n-v(E)  p(x.u.E) 


(19) 


(20) 


n-u  0(x,u,E) 


(21) 
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This  sane  relationship  can  be  obtained  by  considering  the 

number  of  particles  which  cross  the  area  per  unit  time.  For  a 

time  interval  dt,  all  particles  closer  than  v(E)  dt  will  cross 

the  surface  (ignoring  interactions).  The  volume  of  particles 

with  direction  u  that  cross  the  surface  is  a  tube  of  length 

v(E)  dt  and  cross  sectional  area  dA  »  dA  n*u  where  dA  is  the 

n -  n 

differential  area  perpendicular  to  n. 

dV  »  dA  n*u  v(E)  dt  (22) 

n  —  — 


Figure  6.  Flux  Tube  Geometry 


multiplying  the  volume  by  the  particle  density  gives  the  total 
number  of  particles  crossing  the  area  dAn>  Dividing  by  the  area 
and  the  time  interval  dt  yields  the  total  number  of  particles 
crossing  per  unit  area  and  per  unit  time 

*  o(i»JJ.E)  dV/dA^  dt  *  jn  •  u  v{E)  o(x.u.E) 

■  H’H  (23) 

Thus,  the  number  of  particles  crossing  a  surface  with  orienta¬ 
tion  n  is  equal  to  the  flux  only  for  particle  directions  u 
parallel  to  n  so  that  n-u  »  1.  For  all  other  particle  directions, 
the  particles  crossing  per  unit  area  perpendicular  to  n  is  the 
flux  reduced  by  the  factor  n-u. 
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For  space  radiation  problems,  the  initial  population  of 
particles  is  specified  by  the  particle  current  entering  a  surface 
which  surrounds  the  vehicle 


ill! 

4ir 


u*n 


(24) 


The  number  of  particles  entering  this  surface  per  unit  area  is 
giveo  by  an  Integration  over  entering  directions 


Jn(x,E) 


u*n>0 


(25) 


Using  the  cosine  of  the  polar  angle  measured  from  n,  u,  and  the 
azimuthal  angle  measured  around  n.  9, 


jn(x.E) 


2«  1 
l  l 


3-CE>  u 


du  d9  » 

4 


(26) 


Thus,  the  number  of  particles  entering  the  surface  per  unit 
area  is  4/4.  Similarly,  in  the  absence  of  interactions  inside  the 
surface,  the  number  of  particles  leaving  the  surface  per  unit  area 
is  4/4  for  a  net  flow  of  zero.  Both  the  entrance  and  exit  currents 
have  a  cosine  distribution  relative  to  the  surface  normal. 

TRANSPORT  KERNELS 


The  response  at  a  component  is  determined  by  an  integration 
of  the  form 


R 


S°(P)  M(P.P)  D°(P* )  dP  dP* 


(27) 
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where  S°(P)  denotes  the  density  of  source  particles,  M(P,P') 
is  the  flux  at  Pf  per  unit  source  at  P  and  D°(P')  is  the  response 
per  unit  flux  at  P'.  An  abbreviated  notation  has  been  used  for 
particle  coordinates 


P  =  x  »u  ,E 


P'  =  x' ,u ' ,E* 


dP  =  dx  du  dE 


dP '  =  dx'  du '  dE 1 


(28) 


The  discussion  is  now  limited  to  a  fixed  source  with  a 
separable  energy  dependence 

S°(x,u,E)  =  S°(x,u)  f (E)  (29) 

where  f(E)  is  the  normalized  source  spectrum,  and  a  point 
isotropic  detector 


D°( x' ,u ' ,E' )  =  6(x'  -  Xq)  S(E') 


(30) 


where  6  is  the  Dirac  delta  function  and  g(E')  is  the  energy 
dependence  of  the  response  per  unit  flux. 


For  a  uniform  infinite  medium,  the  total  response  can  now 
be  written  as 


R 


f f S°(3.H) 


M(x,u,x  ) 

—  — 


dx  du 


(31) 


where  the  kernal  M  is  the  response  at  per  source  particle 
emitted  at  x  in  the  direction  u  integrated  over  source  energy, 
and  energy  and  direction  at  the  detector 
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1 


MCs.u.Xq)  MCx.u.E.Xq.u'.E’)  g(E')  dE'  du'  dE  (32) 

M  can  be  represented  by 


=  K(r,u-u')/r2  r  *  |x  -  x 


-o  -i 


=  (X  -  X)/3 


(33) 


since  only  the  source  and  detector  separation  r  and  relative 
orientation  u.u'  are  pertinant  parameters  for  an  infinite  medium. 

The  response  kernel  K  can  be  determined  by  any  number  of 
methods,  e.g.,  Monte  Carlo.  For  the  following  discussions,  the 
method  of  calculation  and  the  actual  variation  with  distance  and 
direction  is  not  required. 

The  volume  spatial  integration  over  the  source  for  the 
integral  response  can  be  written  in  spherical  coordinates 


o 

dx  °  r  dr  du' 


(34) 


with  an  origin  at  the  detector  point  x^.  The  R  integral  becomes 
R  “  ///**  •  ru'  ,u)  K(r,u-u')  du'  du  dr  (35) 


For  a  surface  source,  the  source  density  can  be  written 


S°(x,u)  «  S°(x,u)  6(p) 


(36) 
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where  6 ( p )  is  a  one-dimensional  delta  function  in  a  thickness 
variable  p  measured  perpendicular  to  the  surface,  i.e,,  measured 
along  n.  Using  a  delta  function  transformation 

6(p)  -  o(r  -  s)  g  =  «(r  -  s)  (37) 

where  s  is  the  distance  to  the  surface  along  -u’ .  Substituting 
into  the  integral  response  equation  and  performing  the  integra¬ 
tion  over  r  yields 

R  «  ffs0^  -  su'  ,u)  K(s,u/u')  (38) 

If  the  source  per  unit  area  is  spatially  constant 


S"(x  -su' , u )  »  S°(u) 


—o 


•fj*w  «•*•*’> 


(39) 


If  the  source  is  azimuthal ly  symmetric  around  n 


S°(u)  -  3°(n-u)  /  2t? 


R 


1 


JJ S°(n*u) 


K(s,u*u’ ) 

- n  -u 


(40) 


If  the  angular  dependence  relative  to  the  surface  normal  is 

cosine  and  the  total  emission  in  the  forward  hemisphere  is  i 
2 

particles/cro  ,  then 
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S°(n*u)  =  n *u  /  2 

E  =  li  //s-iSts.M’)  ^r-  (41) 

If  particles  transport  straight  ahead^(no  deflections)  then 
K(s,u*u’)  =  K(s)  6(u-u'  -  1)  /  2tt 

*  K(s)  6(u  -  u')  (42) 


^n*u  K(s,u*u’)  du  *  K(s)  n -u' 
'  8  -  4?  /' 


K(s)  du’ 


(43) 


If  the  uniform  infinite  medium  is  a  void 


K(s)  =  1 
R 


si  fAa 


’  -  1 


(44) 


For  general  problems,  K  and  S  are  expanded  as  Legendre 
series 


K(s,u*u' )  -  E  (2i*l)  K,(s)  PH (u-u ' 
*  i  »  0  1  1 - 


) 


(45) 


S°( 


X. 


4*  ^  (2j  +  1)  s”  P,(u-n) 

j  «  0  J  J  — 


(46) 
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!'■  i 


: 


where  P^Cu)  is  the  kth  Legendre  polynomial  and  the  expansion 
coefficients  are 


/ 


K^s)  »  2t r  J  K(s,y )  P^y)  dy 


(47) 


Sj  =  2tt  J  S°(u)  Pj(u-n)  d(u-n) 


(48) 


-1 


Substitution  of  the  series  into  the  response  integral  yields 


E  =  §  S°/Ki(s> /p j«a-D  VH-a)  («) 


The  aadition  theorem  for  Legendre  polynomials’ allows  writing 


PlCu-u’)  55  pi(E*E)  P1(n-u* )  +  f (8) 


(50) 


where  0  is  an  azimuthal  angle  measured  around  u  and  f(0)  is  2ir 
periodic.  Integrating  over  u 


[Pi(u-n)  P1(n-u')  +  f(0)J  du 


4tt 


JTT~J  Pi(n-u')  6 


ij 


(51) 


where  6^  is  the  Krcnecker  delta,  and  the  integral  response  is 


*■  i?o  S?/Ki'=) 


(52) 
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For  a  hypothetical  cosine  surface  source  emitting  i  particle  in 
the  forward  hemisphere,  and  -£  particle  in  the  backward  hemisphere 

S°(u)  =  u-n  /  4tt  Sj  =  6j  1  /  3  (53) 

i.e.,  only  the  first  expansion  coefficient  is  nonzero.  The 
integral  response  becomes 


B  =  4?  /K  1<«> 

*  4?  / *>’ 


(54) 


(55) 


since  P^(n-u')  =  n*u' . 

For  a  point  in  the  center  of  a  sphere,  s  =  t  (the  radius) 
for  all  directions  so  that 


R.(t)  =  Kx(t) 


(56) 


For  a  point  midway  between  two  planes  separated  by  2t,  s  =  t/u 
where  u  is  the  cosine  of  the  angle  measured  from  the  normal  to 
the  planes,  and 


R„(t)  =  2 


1 

4ir 


Kx(t/u)  dy  00 


CO 


ft 


(57) 
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V.  ANOTHER  SECTORING  KERNEL 


The  dose  at  a  point  at  radius  a  inside  a  solid  sphere  of 
radius  b  =  a  +  t,  is  denoted  by  RQ(a,t),  where  the  cosine  source 
is  applied  at  the  larger  radius,  a  +  t,  and  the  sphere  is  a 
subset  of  a  uniform  infinite  medium.  It  follows  from  Equation 


52  that 


Ro(a’t)!34i  f  Ki(s)  du’ 
4ir 


£  f  Kx(s(u))  dp 
ll 


(60) 


where  the  law  of  cosines  yields 


2  2  2 

b  =  a  +  s  +  2  a  sy  —ys  =»  -ya  + 


Vb2  -  -2 


a  (1  -  (61) 


Assuming  the  sphere  is  hollow  for  r  <  a  and  that  the  cavity 
has  little  effect  on  the  inward  dose,  kernel 


1 

Kx(s)  dy  +  £ 

inward  directed  particles  traversing 
particles  cavity 


Rq(  a»t  j 


/ 


/ 


Kx(s)  dy 


(62) 


For  the  particles  traversing  the  cavity,  ■'■he  distance  through 
material  is  now 


-ya  + 

>/b2  -  a2  (1 

V> 

ya  + 

Vb2  -a2  (1  - 

u2) 

-2ya) 


(63) 


so  that  the  upper  and  lower  hemisphere  integrations  are  equal 
and 


X 

RQ(a,t)  •J  K1(s(y))  dy 


(64) 


The  angular  integration  is  transformed  to  an  integration 
over  distance  using 


,  .  r  2  2  2l  .2  2  ,  2  ,2  ,  2 

Z&L  =  zA  b ,  .  *  -  s-.  I  =  b-  -3  ..  .+  .s.,  ,  (65) 

ds  ds  L  2as  J 


2as2  2as2 


where 


h2  -  b2  -  a2  «  (t  +  a)2  -  a2  »  t2  +  2  a  t 


(66) 


so  that 


2  n  r 

R-(a.t)  -  —  f  K,(s)  “I  +  —  /  K,(s)  ds  (67) 
°  2a/1  s2  2a/1 


Assuming 


dK 


KX(S)  *  K.(S)  -  8  gj-  (8) 


(68) 
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still  holds,  integration  by  parts  yields 


\ 


if 


R  <a,t}  «  K=(t)  +  ^  t  K=(t)  -  h  K_(b)  +  /  KJs)  ds  (69) 


+  f  KJs)  dsl 

*4  J 


For  large  a, 


R0(«,t)  =  K=(t) 


(70) 


as  expected;  i.e.,  the  inside  surface  of  the  spherical  shell 
looks  like  a  slab  to  the  transporting  particles. 

For  an  inside  radius  much  less  than  the  shell  thickness, 
a  «  t 

h  =  V?  +  2  a  t  =  t  yjl  +  2  a/t  =  t(l  +  *(2  a/t) ) 


=»  t  +  a 


dK_ 


Kc(s)  =>  Ks(t)  +  i (t)  (s  -  t),  first  term  expansion 


dK= 

K^h)  =  Kja(t)  +  a  -pr-(t)  since  h  -  t  ®  a 


(71) 


Substitution  and  evaluation  of  the  integral  yields 

dKm  a  dK= 

RQ(a,t)  ®  K,(t)  -  t  ^-(t)  -  -g  ^f-(t) 


(72) 


and  for  a  *  0, 


dK. 


Ro(0,t)  *  K„(t)  -  t  -^(t)  -  K.(t) 


(73) 


the  solid  sphere  kernel,  as  expected. 
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The  forward/backward  symmetry  of  the  angular  flux  on  the 
cavity  surface  can  be  seen  from  the  geometry.  For  every  particle 
entering  the  cavity  in  the  direction  u  there  is  an  equivalent 
particle  hitting  the  entrance  point  with  direction  -u. 

The  dose  inside  the  spherical  cavity  can  be  obtained  by  an 
angular  dose  current  integration,  where 

o 

JQ(a,t,u)  =  J0(a,t,-p)  dose/steradian/cm  (74) 


In  fact,  for  large  inner  radii 


J0(°°,t,u)  =  J=(t,u) 


(75) 


where  J=(t,y)  is  the  angular  dose  current  in  the  semi-infinite 
slab  geometry,  with  or  without  a  void  at  the  detector  plane. 


In  terms  of  angular  dose  current 

1 

RQ(a,t)  ■  4ir  f  Jo(a,t,u)  ^ 


(76) 


R„(t) 


4tt 


/ 


J.(t.y) 


lii 


(77) 


since  the  angular  flux  and  current  are  related  by 


J(u)  ■  y$(u) 


or  $(y)  «  J(u)/u 


(78) 
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Assuming  the  transmitted  dose  current  varies  as 


J0(a,t,y' ) 


1_ 

4u 


A(t)  p'a(t) 


(79) 


and  if  RQ(r,a,t)  denotes  the  dose  at  radius  r  in  a  spherical 

cavity  of  radius  a  >  r,  the  dose  can  be  calculated  as 

1 

R0(r,a,t)  =  4ir  f  J0(a,t,uf)  ^  (80) 

*0 


For  r  =  a,  u'  58  u  (Fig.  7c)  and 


R0(a,a,t) 


A(t)  f 

4 


du 


A(t) 

7 


For  r  »  0 ,  u  *  *  1 ,  and 


Ro(0,a,t)  ■  A(t 


A 

)  j  du  “  A( 


t) 


In  the  limit  of  large  a 


RQ(*.®.t)  -  R„(t)  and  Ro(0,«*,t)  -  RQ(t) 

which  implies  that 

A(t)  «  K0(t) 

*U)  -  K0(t)  /  Ka(t) 


(81) 


(82) 


(83) 


(84) 

(85) 


36 


This  transmitted  dose  kernel  is  applied  in  complicated 

geometries  by  assuming  that  materials  traversed  are  sections 

of  a  spherical  shell.  This  requires  dose  attenuation  kernel 

evaluations  based  on  t„  where  t  is  an  estimate  of  the  minimum 

o  o 

path  traversed  by  the  radiation 


tQ  »  t(slant  path)  *  (vt'!!iDSi(1e  ♦  M'lWside5  <86> 


A  combined  kernel  which  accounts  for  both  shell  and  solid 


K(t)  -  K  (t) 


*  ’ 

°s„u> 

a(tQ)  uMto)  “  1 

ar  MO 

1  *•  tJTT  K„(t)  j 

(37) 


shell  shield  term  solid  shield  term 


where  it  is  assuaied  that 


dK.  dK„  KJ  t  ) 


(38) 


t  is  the  distance  through  solid  material  surrounding 
the  dose  point,  zero  if  the  point  is  in  a  cavity 

tQ  is  the  estimated  minimum  distance  through  all  the 
other  shell  like  shields,  aero  if  only  a  so* id  shield 
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(89) 


t  =  t  +  t 
•  o 

a(tQ)  =  K0(t)  /  K=(t0) 

=  1  if  t.  is  zero 

No  attempt  has  been  made  to  include  curvature  effects  for  the 

shell  shields.  For  t  =0.  the  solid  shield  kernel  results  and 

o 

for  t.  =0,  the  angular  dose,  spherical  shell  kernel  accrues. 

The  following  steps  are  used  to  generate  kernel  information 
using  conventional  slab  geometry  codes 

(1)  K_(t)  is  generated  rather  than  K_(t)  to  properly  account 

00 

for  backscattered  particles.  The  problem  has  a  cosine 
source  emitting  J  particle  into  an  infinitely  thick  slab. 

(2)  The  backward  and  forward  angular  doses  (response  weighted 
angular  fluxes)  are  obtained  from  the  same  calculation. 
Assuming  a  Monte  Carlo  calculation 


J5(t,*l<-11°)  •  Ml1-  »„>  S 


backward  traveling 
histories  with  directions 


-lSu  <  -u. 


(90) 


J_(t,u>y_) 


27r(l  -  mo) 


E 


forward  traveling 
histories  with  directions 
uQ  <  u  £  1  (91) 
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(3)  Estimate  the  center  of  spherical  shell  kernel  using 


«J0(t,p>ii0)  *  JQ(t,y<-p0) 


i  J  (t,u>u  )  +  J  (t,p<-u  ) 


KQ(t)  =  4irJo(t,y  =  1) 


r 

J_(t,u>u0)  +  J_(t,U<-U0) 

00  CO  m 


(4)  Calculate  the  ratio  of  spherical  shell  to  slab  geometry 
dose 


o(t)  =  K  (t)  /  KJt)  >  1 


VI.  SECTORING  KERNEL  SELECTION 


Charged  particle  dose  predictions  for  satellite  and  deep 
space  probes  are  usually  performed  by  sectoring  calculations. 

The  accuracy  of  these  predictions  is  controlled  partially  by  one¬ 
dimensional  dose  attenuation  kernels.  This  section  discusses 
kernel  selection  in  the  energy  domains  where  particles  have 
straight  line  paths  and  where  the  particles  have  significant 
angular  deflections. 

Sectoring  calculations  assume  an  explicit  three-dimensional 
analytic  geometry  mockup  of  the  spacecraft.  The  dose,  D,  at 
critical  points  in  the  spacecraft  is  computed  as 

d2  &  (95) 

where  u  is  a  direction  from  the  point  to  the  outside  of  the 

vehicle,  t(u)  is  the  material  thickness  encountered  along  u, 
o 

d  u  is  the  differential  solid  angle,  and  K(t)  is  a  dose  attenua¬ 
tion  kernel. 

This  integration  is  performed  numerically  by  defining  many 
solid  angle  sectors;  e.g.,  I  total  sectors  each  of  size  4tt/I 
steradiaas,  for  each  sector  i  defining  an  average  direction 
u^,  and  evaluating 
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sector 

solid 

angle 


(96) 


i=l 


4tt 


Sum 

over 

sectors 


dose  per 
steradian 

for  thickness  t(u^) 


The  thickness  t(u^)  is  determined  by  a  ray  trace  from  the 
dose  point  along  the  direction  u^.  Three  sources  of  error  in 
this  type  of  analysis  are: 


(1)  The  fineness  of  the  angular  integration  which  determines 
how  well  thin  vehicle  sections  are  resolved  (or  if 

they  are  seen  at  ail). 

(2)  The  degree  to  which  the  radiation  arriving  from  the 
direction  u.  actually  traverses  the  materials  comprising 
the  thickness  t(ui).  This  approximation  neglects  the 
scatter  and  change  of  direction  as  the  particles  traverse 
the  materials. 

(3)  The  extent  that  simple  geometry  used  in  deriving  the 
attenuation  kernel  K(t)  approximates  the  actual 
geometry  of  the  sector. 


The  first  error  source  will  not  be  discussed  further  here. 
It  can  be  investigated  in  most  codes  by  parameters  controlling 
the  size  of  the  solid  angle  sectors. 

The  latter  two  error  sources  are  related  since  they  both 
require  knowledge  of  how  particles  really  get  from  the  outside 
of  the  vehicle  to  the  dose  point.  A  definitive  answer  to  the 
error  magnitude  is  outside  the  realm  of  sectoring  analysis  and 


requires  numerical  nethods  such  as  adjoint  Monte  Carlo  (which 
provides  a  realistic  simulation  of  all  the  processes  working  on 
the  particles  as  they  progress  towards  the  detector  point). 

However,  it  is  possible  by  examining  some  typical  one-dimen¬ 
sional  geometry  problems,  to  get  some  insight  into  the  potential 
magnitude  of  the  errors  in  sectoring  analysis  and  to  select  one 
dimensional  kernels  that  will  provide  both  upper  and  lower 
bounds  on  the  dose  for  most  three-dimensional  geometries. 

The  capability  of  bracketing  the  true  dose — by  sectoring 
with  upper  and  lower  bound  kernels — will  alleviate  the  accuracy 
question  for  many  dose  points.  When  the  upper/lower  bound 
kernels  are  in  substantial  agreement,  differing  by  less  than  a 
factor  of  2,  then  a  good  dose  prediction  can  be  assumed  (using 
the  upper  bound  kernel  for  conservatism  in  design).  When  there 
is  substantial  disagreement  between  the  upper  and  lower  bound 
dose  estimates — more  than  a  factor  of  2  disagreement — further 
analysis  is  required. 

Additional  analysis  can  include  adjoint  Monte  Carlo,  or  an 
in-depth  review  of  the  dose  point  and  its  surroundings.  By 
reviewing  the  local  geometry  of  the  dose  point,  it  may  be  possible 
to  identify  which  of  the  kernel  estimates  is  expected  to  best 
represent  a  particular  dose  point. 

The  remainder  of  this  document  covers  one-dimensional  kernels 
in  some  depth  with  the  following  intent: 
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(1)  Identify  the  correct  kernels  for  specific  particle  types; 

(2)  identify  the  difference  in  dose  predicted  by  different 
kernels  and  the  source  of  that  difference,  and 

(3)  provide  some  arguments  for  using  specific  kernels  in 
particular  circumstances. 

STRAIGHT  AHEAD  MODELS 

Sectoring  is  a  straight  ahead  transport  model  since  the 
particle  transport  is  characterized  only  by  the  materials — t(u) — 
encountered  along  a  specific  direction  u. 

The  straight  ahead  method  is  reliable  for  some  analyses, 

e.g. , 

(1)  Dose  from  cosmic  rays  and  protons,  and 

(2)  dose  from  electrons  with  energies  exceeding  100  MeV. 

Some  sectoring  codes  do  a  complete  particle  energy  analysis 
in  the  straight  ahead  approximation.  The  integral  evaluated  is, 
for  a  given  sector 

00  00 

K(t)  »  f  f  $0(E')  M(E'-E.t)  Dq(E)  dE 1  dE  (97) 

O  *o 

where  K(t)  is  the  dose  at  thickness  t,  $o(E')  is  the  external 

2 

environment  (particles/cm  •  MeV),  M(E t)  is  the  distribution 
in  energy  E  at  the  detector  point  for  a  particle  that  starts  with 
energy  E'  and  traverses  material  thickness  t,  and  DQ(E)  is  the 
detector  response  (dose)  to  a  particle  of  energy  E, 
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!^r.«Niwn 


J 

1 

t 


In  the  straight  ahead  approximation,  particles  lose  energy 
at  a  rate  per  unit  path  length  s  given  by  L(E),  i.e. 


dE 

ds 


(98) 


Therefore,  after  traversing  the  thickness  t,  a  particle 
that  starts  at  energy  E'  at  the  exterior  of  the  satellite  reaches 
the  detector  point  with  energy  E  given  by 


E 


f ( E ’ , t )  =  E’ 


L(E(s) )  ds 


(99) 


Neglecting  all  other  processes, the  basic  transport  kernel 
M(E  -*•  E,t)  becomes 


M(E'-*>  E.t)  »  6(E  -  f(E',t))  (100) 

where  <5  is  the  Dirac  delta  function,  so  that  the  dose  K(t) 
becomes 

00  00 

K(t)  =  //  d>0(E ’ )  6 ( E  -  f(E\t))  D0(E)  dEdE' 

"o  o 

CO 

- J  $Q(E ' )  Do(f(E’,t))  dE’  (101) 

o 

The  simplest  and  most  often  used  method  of  evaluating  this 
integral  is  to  subdivide  the  energy  domain  into  intervals  called 
groups  and  evaluate 
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(102) 


K(t,  .  E 

g 

groups 


*0CEg)AEg 

particles  in 
group  g 
impinging  on 
vehicle 


VEg(t)) 

response  to  these 
particles  at  their 
degraded  energy  at 
the  detector 


where  g  is  a  group  subscript,  E’  is  the  average  energy  of  the 

o 

particles  in  group  g,  AEg  is  the  group  width,  Eg(t)  is  the  energy 
the  particles  degrade  to  after  traversing  the  thickness  t 


Vt) 


L(E(s))ds 


(103) 


and  DQ(Eg(t))  is  the  response  to  these  particles. 

This  simple  method  has  one  major  problem.  Since  charged 
particles  range  out,  all  groups  above  some  critical  energy  make 
it  to  the  detector  point  while  all  those  below  that  energy  stop 
before  getting  there.  This  yes/no  on  groups  can  cause  bumpy 
dose  predictions  unless  a  fine  group  structure  is  used. 


An  alternate  way  of  evaluating  the  straight  ahead  dose  is 
presented  here  because  it  eliminates  the  bump  problem  and  is  a 
major  item  that  distinguishes  adjoint  Monte  Carlo  from  conven¬ 
tional  particle  simulation  Monte  Carlo. 


In  the  integral 


K(t)  = 


♦o<E’)  D0(f(E\t))  dE* 


(104) 
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make  the  change  of  variables 


E  ■  f(.E',t)4-4E’  =  f-1(E  ,t) 


(105) 


i.e. 


E'  =  f "1(E , t) 


E  + 


L(E ' ( s) )  is 


(106) 


then 


K(t) 


(E , t) )  Dq(E) 


dE 


(107) 


where jdE' /dE|  is  the  Jacobian  of  the  variable  transformation.  A 
careful  evaluation  of  the  Jacobian  yields,  for  a  single  material, 


dE’ 

dE 


(108) 


Therefore,  the  dose  is  also  given  by  an  integration  over  energy 
at  the  detector  point  in  the  form 


K(t)  = 


00 

f  $0(f'1(E’t)) 

O 


Dq(E)  dE 


(109) 


which  can  be  evaluated  using  group  terminology  as 


K(t) 


E 

s 

groups 


L(  E '  ( t ) ) 

W*»  -rf=—  Do<V“g 


(110) 


Now  the  groups  are  defined  at  the  detector,  E  is  the  average 

& 

energy  of  the  particle  that  arrive  at  the  detector  (preassigned) 


in  group  g,  E'(t)  is  the  energy  these  particles  need  to  hit  the 
6 

satellite  with  to  arrive  at  the  detector  with  energy  E  ,  i.e., 

& 

t 

^g(t)  »  !g  +  f  L(E '  (s) )  ds  (111) 

*o 

This  scheme  eliminates  any  on/off  bumps  at  the  detector. 

The  major  difference  in  this  scheme  is  that  particles  speed 
up  in  backtracking  from  the  detector  point  to  the  source,  and 
that  a  factor  L(E')/L(E)  is  introduced  into  the  equation  for 
dose.  The  L(E')/L(E)  factor  accounts  for  the  contraction/ 
expansion  of  particle  group  width  during  flight  through  the 
vehicle  materials. 

Generally,  stopping  powers  L(E)  get  larger  at  lower  energies. 
This  causes  group  widths  to  expand  as  particle  travel  from  the 
source  (outside  the  satellite)  to  the  detector.  Conversely,  in 
tracking  particle  groups  from  the  detector  to  the  source,  groups 
get  narrower. 

The  factor  L(E ' ( t ) )/L(E)  is  the  ratio  of  group  width  at  the 
source  to  group  width  at  the  detector  point  for  an  infinitesimal 
width  group.  (Multiple  materials  require  products  of  L(E')/L(E) 
terms  corresponding  to  the  entrance  energy  and  exit  energy  for 
each  material  layer.) 

Most  sectoring  codes  do  not  evaluate  transport  integrals 
explicitly.  The  integration  is  performed  in  a  one-dimensional 
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code  and  results  are  supplied  to  the  sectoring  code  as  tabular 
data  t^,  K(t^),  i  *  1,  2,  ...  .  The  dose  for  a  particular  sector 
thickness  K(t(u))  is  then  obtained  by  interpolation  on  thickness 
t. 

As  mentioned  previously,  the  straight  ahead  method  is 
accurate  for  heavy  charged  particles  and  very  high  energy  elec¬ 
trons  (>  100  MeV).  However,  some  care  must  be  used  in  selecting 
the  kernel  K(t)  since  there  are  a  variety  of  one-dimensional  codes 
that  can  generate  kernels  for  different  source  and  material  geome¬ 
tries.  The  applicable  kernel  is  one  of  the  following  (they  are 
all  equal  numerically  because  of  the  straight  ahead  model): 

(1)  Dose  at  depth  t  in  a  slab  for  a  source  <t>Q(E)  normally 
incident  on  the  slab, 

(2)  dose  at  the  center  of  a  solid  sphere  of  radius  t  for 
any  source  that  has  an  angular  strength  of  4>0(E)/4tt 
particles  per  unit  area  and  per  unit  steradian  normal 
to  the  sphere  surface  (particles  incident  at  angles 
other  than  normal  never  get  to  the  sphere  center), 

(3)  dose  at  the  center  of  a  spherical  shell  of  thickness  t 
with  the  same  source  conditions  as  2,  and 

2 

(4)  4wt  tiroes  the  dose  at  the  outside  surface  of  a  sphere 
of  outer  radius  t  with  a  point  source  of  $0(E)  particles 
MeV  at  the  center  of  the  sphere. 


48 


One  kernel  that  must  not  be  used  is  the  dose  at  the  center 
of  a  slab  of  thickness  2t  where  that  slab  is  embedded  in  an 
isotropic  flux  field  of  <{>  (E). 

Denoting  this  kernel  by  K_(t)  and  the  correct  kernel  by 
K.(t)  (the  subscript  look  like  a  slab  =  and  a  sphere  . ),  it  is 
seen  that 


K=(t) 


f  K.(t(u)) 
J  4tt 

4ir 


(112) 


2 

Writing  d  u  as  da  dy  were  a  is  azimuth  around  the  slab 
normal  and  y  is  the  cosine  of  the  angle  between  u  and  the  slab 
normal  then 


Ka(t) 


L  L  da  du 

i  ^  i 

2,  .  2j  du  . J  K.(t/u)  du 


(113) 


Since  dose  usually  decreases  with  thickness,  K=(t)  is  less 
than  K  (t).  Therefore,  the  use  of  a  slab  kernel  in  a  sectoring 
code  will  result  in  an  underestimate  of  the  true  dose.  (Imagine 
applying  the  sectoring  code  to  the  sphere  or  slab — neither  geometry 
gets  the  right  answer  if  Ka(t)  is  the  input  kernel.) 
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Making  the  variable  transformation 


x  -  t/y  — ^  y  *  t/x  and  dy  =  t-t/x")  dv 


K  (t)  »  t 


00 

/ 


(*)  ^ 


(114) 


(115) 


Differentiating  with  respect  to  shield  thickness  t 


1 

j 

00 

at  U‘)  -  /  m*>  %  ♦ 1 

i 

w 

to 

X 

j 

J%  x  1 

x 

x=t 

j 

-  1  M‘>  -  ^ 

(116) 

1 

thus 

s.(t)  =  K.(t)  -  t  K.(t) 

(117) 

A  comment  on  angular  distributions  is  warranted  because 
it  leads  naturally  to  the  next  subject.  The  slab  kernel  can 
be  written  as 

K  ft)  -  /  K_(t,y)  dy 


/-• 


(118} 


where  Ks(t,u)  is  the  dose  per  steradian  for  particles  arriving 
at  the  detector  at  an  angle  cos”*u  relative  to  the  slab  normal. 
Since 


K_(t) 


/  «• 


(t/u)  ,*.y 
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(119) 


it  follows  that 


4*K=(t,y)  -  K.(t/p)  (120) 

or 

K,(t,p)  =  K.(t/w)  (121) 

Therefore,  for  straight  ahead  dose  attenuation  kernels 
K.(t)  that  decrease  with  thickness,  the  angular  distribution  o * 
Km(t,u)  gets  sharply  peaked  along  the  slab  normal.  (K#(t/y)<<K,(t) 
for  values  of  y  less  than  !•)  Said  another  way,  only  those 
particles  that  are  incident  on  the  slab  in  some  area  around  the 
normal  contribute  to  the  dose. 

ELECTRONS 

The  transport  of  electrons  below  10  MeV  (typical  earth 
environments)  is  not  described  by  straight  ahead  models.  The 
transmission  of  these  electrons  through  materials  is  very  close 
to  a  diffusion  process;  eg.,  the  angular  distribution  of 
particle  transmitted  through  surfaces  is  almost  isotropie. 

Moreover,  the  relative  level  of  dose  values  is  determined  by 
the  actual  thickness  of  components  such  as  box  walls  rather 
than  the  slant  thickness  seen  by  a  sectoring  code. 

One  indication  of  the  effect  of  this  diffuse  transport  can 
be  seen  in  simple  spherical  geometries.  The  inside  surface  of 
a  spherical  shell  of  large  radius  looks  like  a  slab  to  local 
incoming  electrons  (the  electron  range  in  the  shell  material  is 
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much  less  than  the  shell  radius).  If  the  transmitted  electrons 
retain  their  isotropic  distribution,  the  dose  at  the  center  of 
the  sphere,  Ko(t),  and  everywhere  inside  the  shell  is  constant 
and  euqal  to  the  slab  dose  K_(t). 

In  practice,  the  electrons  transmitted  through  a  slab  of 
thickness,  t,  do  become  forward  peaked,  but  not  severely.  The 
angular  transmission  can  be  approximated  by 

K=(t,y)  =  4irKo(t)ya(t)-1  (122) 


where 

a(t)  =  Ko(t)/K=(t) 


(123) 


and  is  typically  in  the  range  1  (isotropic)  to  1.5.  (Other 
angular  distributions  can  be  used  but  require  more  parameters. 

The  above  distribution  suffers,  foe  instance,  by  predicting 
K_(t,u=0)  =  0  if  a(t)  >  I).  Thus,  in  reality,  the  electron 
dose  inside  the  sphere  varies  from  X_(t)  at  the  inside  surface 
to  K^(t)  at  the  center,  usually  less  than  50  percent. 

This  small  variation  of  dose  inside  a  spherical  shell  for 
electrons  should  be  contrasted  with  the  dose  variation  obtained 
for  protons  using  straight  ahead  models.  Since 

K.(t)  =  K„(t)  -  t  ^  K=(t)  (134) 

the  variation  from  inside  surface  (K_)  to  center  (Ko)  is  often 
an  order  of  magnitude  or  more  for  protons. 
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Dose  enhancement  is  also  seen  for  electrons  as  the  inner 
radius  of  the  spherical  shell  is  shrunk  to  values  on  the  order  of 
the  electron  range  or  less— the  dose  remains  roughly  constant 
inside  the  cavity  hut  is  higher  for  the  smaller  cavity  radius. 
Specifically,  it  has  been  shown  that  for  a  detector  point  inside 
a  solid  spherical  shield,  that  the  electron  dose  is  given  by 
K.(t)  and  not  K^Ct)  (remember  that  KQ  is  the  dose  at  the  center 
of  a  very  lar^e, nearly  infinite,  radius  sphere)  where  K(t)  is 
computed  from  K_(t)  -  t  K=(t).  This  result  has  also  been 
demonstrated  by  adjoint  Monte  Carlo  calculations. 

Note  thau  Ko(t)  and  K  (t)  are  not  equal  for  electrons  ( they 
are  for  straight  ahead  models).  In  fact,  since  KQ(t)  is  on  the 
order  of  K=(t)  (or  as  much  as  50  percent  larger),  and  since  K  (t) 
can  be  an  order  of  magnitude  largei  than  K_(t),  the  difference 
between  KQ(t)  and  K(t)  can  also  be  an  order  of  magnitude  for 
electrons. 

Reiterating  the  above,  it  is  seen  that  the  electron  dose 
anywhere  inside  a  large  spherical  cavity  with  shell  thickness 
t  is  approximately  K=(t).  A  sectoring  code  will  yield  this  same 
result  for  the  shell  if  Ki(t>  is  supplied  as  input  and  the  dose 
point  is  located  at  the  center  of  the  sphere.  However,  if  the 
dose  point  is  located  near  the  inner  surface  of  the  shell,  the 
sectoring  code  will  underpredict  the  dose  since  an  evaluation  of 

1 

Ka(t/u)  du  <  K=( t)  (125) 
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is  effectively  obtained.  This  type  of  conservatism  (underpre¬ 
diction)  cannot  be  tolerated  in  design  calculations. 

This  tendency  to  underpredict  cavity  dose  has  been  eliminated 
in  the  SIGMA-II  sectoring  code  by  forcing  the  program  to  estimate 
the  minimum  thickness,  t,  seen  by  electrons.  The  usual  slant 
thickness,  t,  calculated  by  the  sectoring  code  is  multiplied  by 
cosine  of  the  angle  between  the  direction  u  and  the  normal  at 
entrance  and  exits  of  u  in  the  material 


The  use  of  the  minimum  path,  t,  as  the  thickness  parameter 
is  a  much  better  characterization  of  the  true  transport  character¬ 
istic  of  electrons.  The  additional  use  of  the  functional  form 
^a(t)-l  tQ  characterize  departure  from  isotropic  transmission 
through  material  layers  results  in  a  sectoring  formalism  that 
will  reproduce  the  correct  dose  in  two  simple  geometries: 

(1)  dose  inside  a  spherical  shell,  and 

(2)  dose  behind  a  slab  shield. 

However,  the  above  kernel  does  not  account  for  the  dose 
enhancement  seen  in  small  solids  containing  dose  points  (or  in 
small  inner  radius  cavities).  If  the  small  solid  is  viewed 
separ »tely--ignoring  the  effects  of  plates  and  voids  prior  to 
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reaching  the  solid — it  has  been  shown*  that  the  dose  at  points 
inside  the  solid  is  given  by  using  K.(t)  as  the  sectoring  kernel. 

Therefore,  if  K,(t)  is  used  as  a  sectoring  kernel  with  the 
usual  slant  path  thickness, t,  it  will  also  predict  the  correct 
dose  in  two  simple  geometries: 

(1)  dose  inside  a  solid  sphere,  and 

(2)  dose  behind  a  slab  shield. 

The  question  remains  as  to  which  of  these  kernels  is  most 
appropriate  in  design  calculations.  In  the  limiting  extreme 
of  a  spherical  shell  shield  with  a  very  small  inner  radius,  the 
minimum  path  kernel  underpredicts  dose  while  in  the  other 
extreme  of  a  very  large  inner  radius  the  solid  shield  kernel  over 
predicts  dose. 

Since  underestimates  of  dose  can  be  critical  to  a  mission, 
it  would  seem  that  using  the  solid  sphere  kernel  is  best.  How¬ 
ever,  shield  penalties  associated  with  over  conservative  dose 
estimates  can  also  jeopardize  a  mission.  Consistency  with 
earlier  arguments  would  dictate  using  the  solid  sphere  kernel 
and  slant  path  sectoring.  However,  solid  shield  geometries  are 
seldom  encountered  in  practice — each  sector  usually  involves 

*  "shown"  by  explicit  adjoint  Monte  Carlo  calculation  and  by  an 
analytic  derivation  for  an  idealized  problem  that  has,  to 
first  order,  all  the  characteristics  of  the  real  problem  and, 
to  second  order,  cancelling  terms  that  characterize  the  depar¬ 
ture  of  the  idealized  problem  from  the  real  problem. 


material  layers  separated  by  voids.  Calculations  on  several 
typical  points  using  adjoint  Monte  Carlo  correlate  much  better 
with  the  minimum  path  kernel. 

There  is  a  way  to  combine  the  two  kernels  in  a  systematic 
method  that  accounts  for  peaking  in  small  radius  cavities.  This 

combined  kernel  does  not  exist  in  any  present  code  but  could  be 
added  relatively  easily  to  SIGMA-II  for  instance.  The  kernel 
is  discussed  at  the  end  of  this  section. 

The  correct  kernels  for  electron  dose  calculations  can  be 
obtained  from  one-dimensional  Monte  Carlo  codes.  Adjoint  Monte 
Carlo  can  calculate  the  slab,  solid  sphere,  and  spherical  shell 
configurations  exactly. 

Conventional  particle  simulation  Monte  Carlo  is  limited  to 
calculations  in  slab  geometry  (it  is  hard  to  get  particles  to 
pass  through  a  small  volume  at  the  center  of  the  spheres). 
However,  both  the  spherical  shell  and  solid  sphere  kernels  can 
be  obtained  from  these  slab  geometries  codes. 

Specifically,  the  solid  sphere  electron  dose  kernel  can  be 
obtained  by  the  formula 

•  dK_( t) 

K.(t)  -  K.(t)  -  t  -- -~T—  (127) 

The  validity  of  the  formula  has  been  demonstrated  by  explicit 
adjoint  Monte  Carlo  calculations.  (The  formula  is  not  exact, 
but  any  discrepancies  are  comparable  to  the  Monte  Carlo  error 
bars).  The  spherical  shell  dose  kernel  can  also  be  obtained 


from  a  slab  Monte  Carlo  calculation  by  requesting  the  angular 
dose  (dose  per  steradian)  in  the  direction  along  the  slab  normal. 
The  4ir  times  this  angular  dose  is  the  same  as  the  spherical  shell 
dose  (the  shell  looks  locally  like  a  slab  and  the  angular  flux 
normal  to  the  shell  is  the  only  flux  that  reaches  the  center  of 
the  sphere) . 

The  source  mockup  is  critical  to  the  above  calculations. 

The  dose  kernels  are  supposed  to  be  representative  of  a  shield 
embedded  in  an  isotropic  flux  field.  This  isotropic  flux  field 
is  simulated  by  a  cosine  source  on  the  surface  of  the  shield 
with  a  total  particle  source  of  £  for  a  unit  isotropic  flux. 

(It  should  be  noted  that  data  supplied  by  the  National  Space 
Sciences  Data  Center  on  dose  attenuation  is  iK_(t)  i.e.,  for 
i  of  a  particle — cosine  incidence — a  one  sided  slab  result). 

Slab  geometry  Monte  Carlo  calculations  can  be  performed  for 
finite  thickness  slabs  or  infinitely  backed  slabs.  The  better 
kernel  is  assumed  to  be  the  infinitely  backed  kernel  since  it 
includes  a  backscattering  effect  which  also  occurs  in  the 
satellite,  and  since  it  is  conservatively  higher. 

HIGH  ENERGY  ELECTRONS 

The  final  topic  is  a  specific  to  the  Jovian  electron 
environment.  This  environment  includes  electrons  >  100  MeV. 
Therefore,  the  basic  transport  mechanisms  include  combinations 
of  straight  ahead  and  diffusion  mechanisms.  Specifically, 
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electrons  of  high  energy  transport  in  an  almost  straight  ahead 
manner  and  then  diffuse  after  reduction  to  lower  energies. 

One  way  of  handling  this  type  problem  is  to  use  multiple 

kernels  to  describe  the  complete  range  of  energies.  Each  band 

g  would  have  a  thickness  t  specified  such  that  for  thickness,  t, 

© 

less  than  t  straight  ahead  models  are  employed  while  for  t  >  tg 

the  thickness  t'  =  (t  -  t  )  (closest  to  the  detector)  would  use 

© 

minimum  thickness  models.  The  thickness,  t  ,  would  be  defined  by 
comparing  one-dimensional  calculations  with  straight  ahead  and  with 
explicit  transport  models  and  defining  t  as  the  thickness  where 
the  calcualtions  depart  significantly  in  value.  This  logic  does 
not  exist  in  any  sectoring  code  at  this  time. 

A  KERNEL  WITH  CURVATURE  EFFECTS 

The  difference  in  dose  obtained  behind  a  solid  sphere  and 
a  slab  shield  is  geometrically  obvious  in  the  straight  ahead 
approximation. 

For  the  solid  sphere, the  point  at  the  center  of  the  sphere 
gets  particle  arriving  from  the  entire  4it  steradians  and  each 
particle  has  traversed  the  same  thickness,  t.  For  the  slab,  all 
particles  except  those  traversing  a  path  along  the  normal,  tra¬ 
verse  a  path  larger  than  t.  The  net  affect  is  that  only  a  part 
of  the  4ir  steradians  contributes  to  the  dose  at  depth  t. 

The  came  geometric  effect  occurs  for  electrons  even  though 
the  electron  transport  is  more  diffuse.  The  diffuse  nature  of 
the  transport  makes  it  more  difficult— for  the  slab  geometry— 
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to  define  a  precise  cone  boundary  within  which  particles  traverse 
a  thickness  that  allows  them  to  reach  the  detector.  However, 
this  cone  exists  even  though  its  boundaries  are  fuzzy. 

The  effect  of  curvature  on  cavity  dose  is  predicted  using 
the  fact  that  electron  dose  is  fairly  uniform  inside  cavities 
and  that  the  electron  dose  at  the  inner  surface  of  a  cavity  can 
be  predicted  from  K,(t)  using  slant  thicknesses.  Therefore,  it 
is  possible  to  compute  a  correction  factor  that  accounts  for 
the  relative  size  of  cavities  sncountered  during  sectoring 
since 

K.  (t)  *  Kjt)  -  t  K=(t)  (128) 


in  solid  materials  even  when  the  straight  ahead  model  does  not 
apply. 


Assume  a  is  the  inner  radius  of  cavity  and  x  is  the  spheri¬ 
cal  shell  thickness.  Then  the  dose  at  the  inner  surface  of  the 
cavity  is 


K(x)  *  Ka(x)  +  ^  tK^t)  -  hK  ( b )  +  /  K  (x)dx  (129) 


where  h  is  the  maximum  path  in  the  shell 


h 


(130) 


If  a  is  large,  the  evaluation  gives  Ko(t)  and  if  a  is  small, 
the  value  becomes  K  (t).  (Note  that  in  the  above  equation  a, 
t,  x,  and  h  are  all  in  equivalent  units.  Since  a  is  an  actual 
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dimension  of  a  cavity,  material  thicknesses  must  be  expressed  in 
the  same  units  as  a  (cm)  except  when  interpolating  the  kernels). 

Thus,  every  time  a  cavity  is  encountered,  an  enhancement 
factor  can  be  computed  for  that  cavity  as 


enhancement 

factor 


K_(>r) 

kZTTT 


K,(>h) 

■  u? r 


(131) 


where 


K=(>x) 


a 

l 


K=(y)  dx; 


(132) 


Thus,  an  equation  exists  to  predict  the  dose  enhancement 
from  solid  shields  surrounding  dose  points.  This  equation  holds 
specifically  for  spherical  geometry  surrounding  the  detector 
point.  A  critical  parameter  is  the  radius  of  the  cavity. 

This  section  outlines  a  scheme  for  calculating  an  effective 
radius  of  cavities.  No  claim  for  correctness  of  the  method  is 
made  except  that  it  yields  the  correct  results  for  the  limiting 
extremes  of  slab  and  spherical  shields. 


The  curvature  in  a  plane  containing  Xj  and  x^  - 
either  (x,y),  or  (y,z),  or  (z,x)  —  is  given  by 


-  x^ ,  Xj  denotes 


K 


ij 


dV 


<Lx‘ 


1  + 


(?) 


3/2 


(133) 
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Therefore,  when  a  ray  crosses  a  surface,  compute  the  curva¬ 
ture  in  the  xy,  yz,  and  zx  planes  from  the  surface  equation  and 
define  an  effective  curvature  as  the  average  of  these  three 
values  (or  the  square  root  of  the  average  squares  ...)•  Define 
the  cavity  radius  a,  as  the  reciprocal  of  this  average  curvature. 

This  process  yields  the  following  results  for  specific 
surfaces : 

plane  surface  K  =  0  a  =  <»,  correct 

spherical  surface  K  =  1/a  a  =  a,  correct 

cylindrical  surface  K  =  l/2a  a  »  2a,  ? 

Note  that  the  correct  curvature  is  obtained  for  slab  and 
spherical  geometries.  Therefore,  this  prescription  for  cavity 
radius,  in  conjunction  with  the  equation  for  dose  enhancement 
as  the  cavity  radius  decreases,  will  give  a  kernel  prescription 
that  is  correct  for  slabs,  spherical  she1 Is,  and  solid  spheres. 
The  fact  that  the  kernel  will  correctly  predict  these  three 
geometries,  including  cavity  radius  effects,  is  a  strong  argue- 
ment  for  its  use  in  future  calculations. 
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APPENDIX  A 


MONOENERGETIC  ELECTRON  ATTENUATION  DATA 

Dose  attenuation  curves  are  presented  for  monoenergetic 
electrons  and  their  secondary  bremsstrahlung.  Electron  energies 
are  8,  5,  3,  2,  1,  0.5,  0.25,  and  0.1  MeV.  Shield  materials  are 
hydrogen,  helium,  beryllium,  aluminum,  copper,  tungsten,  and 
uranium. 

Attenuation  curves  are  presented  for  three  shield  geometries, 

n 

dose  between  two  slabs  each  of  thickness  t,D_(t),  dose  at  the  \ 
center  of  a  spherical  shell  of  very  large  inner  radius  DQ(t), 
and  dose  at  the  center  of  a  solid  spherical  shield  D  (t).  Each 
of  the  shields  is  imbedded  in  a  unit  isotropic  free  space  fluence. 

Electron  and  bremsstrahlung  transport  were  simulated  by  a 
combination  of  analytic  and  Monte  Carlo  methods.  Fluctuations 
in  the  computer  results  were  not  smoothed  prior  to  plotting. 
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